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A hierarchy of length scales for weak solutions of the 
three-dimensional Navier-Stokes equations 

J. D. Gibborfl 
Abstract 

Moments of the vorticity are used to define and estimate a hier- 
archy of time-averaged inverse length scales for weak solutions of the 
three-dimensional, incompressible Navier-Stokes equations on a peri- 
odic box. The estimate for the smallest of these inverse scales coincides 
with the inverse Kolmogorov length but thereafter the exponents of 
the Reynolds number rise rapidly for correspondingly higher moments. 
The implications of these results for the computational resolution of 
small scale vortical structures are discussed. 
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1 Introduction 

Resolution issues in computations of solutions of the three-dimensional 
Navier-Stokes equations are not only closely associated with the problem 
of regularity but they also raise the question of how resolution length scales 
can be defined and estimated. The Kolmogorov school of statistical turbu- 
lence suggests that a system of volume has a cut-off in its —5/3 energy 
spectrum at kc = A^^ ^ L~^Re^^'^, which is known as the inverse Kolmogorov 
length. The wave-numbers k>kc are considered to lie in what is called the 
dissipation range [UlSllS]. Significant energy lying in this range can provoke 
intermittent events in the vorticity and strain fields characterized by violent, 
spiky departures away from space-time averages whose corresponding statis- 
tics appear to be non-Gaussian in character [H El El El |8] , although intermit- 
tent events may also be associated with the inertial range (k < kc). Whether 
significant energy actually cascades down to the micro/nano-scales where 
the equation fails to be a valid model is intimately entwined not only with 
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the open question of regularity but also with the role of the Navier-Stokes 
equations as a limit of kinetic theory [9l [10] . This phenomenon continues to 
pose severe computational challenges [HI [12l [13l [ll] . In statistical physics 
the objects that are used to study intermittency are the ensemble-averaged 
velocity structure functions 

(|^(a; + r)-n(a;)|^'),„~/^ (1.1) 

the departure of whose exponents Cp from lineaiH is thought to be caused by 
inertial range intermittent behaviour [Tl|2l[3]. It is clear, however, that these 
structure functions are not best suited for Navier-Stokes analysis : the task of 
this paper is to discuss what could replace these in the Navier-Stokes context 
and what information could be gleaned from them. While higher gradients 
of the velocity field would undoubtedly capture intermittent behaviour, they 
would be unreachable computationally for all practical purposes. A better 
diagnostic of spikiness in Navier-Stokes solutions is a sequence of L^-norms, 
or higher moments, of the vorticity a; = curlti, defined through the set of 
frequencies {p = 2m for m> 1) 

nm{t)=iL-' J^\u;\^"'dv\ , (1.2) 

where the domain V=[0,L]^ is taken to be periodic. The basic frequency 
associated with the domain is given by ijjq = uL~^ . ^\ is the enstrophy 
per unit volume which is related to the energy dissipation rate, whereas the 
higher moments will naturally pick up events at smaller scales. 

The setting is the incompressible (divw = 0), forced, three-dimensional 
Navier-Stokes equations for the velocity field u{x, t) 

dtU + u-Vu = iyAu-Vp + f{x). (1.3) 

Traditionally, most estimates in Navier-Stokes analysis have been found in 
terms of the Grashof number Gr, which is expressed in terms of the root 
mean square (/rms = -^"^ll/lli) °f divergence-free forcing f{x) (see [151 
[T6l [T71 [T8] ) but it would be more helpful to express these in terms of the 
Reynolds number Re to facilitate comparison with the results of statistical 
physics. The definitions of Gr and Re are 

Gr = L^fr„,si^~^, Re = UQLu-\ (1.4) 

^Kolmogorov predicted a linear relation between C,p and p : departure from this is called 
'anomalous scaling' and is usually manifest by C,p lying on a concave curve below linear 
for p>3 [T|. The two coincide for p = 3. 
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Doering and Foias [19] used the idea of defining Uq as 

U^ = L-W\u\\l)^ (1.5) 
where the time average ( • )^ over an interval [0, T] is defined by 

1 

(s'(-))t = 1™sup- / g{T)dT. (1.6) 
9(0) ^ Jo 

Clearly, Gr is fixed provided / is L^-bounded, while Re is the system re- 
sponse to this forcing. A brief look at Leray's energy inequality shows why 
this definition of Uq is of value 



1^ 



u\^dV<-v \u\^dV+\\f\\2\\u\\2, (1.7) 



because it leads to 

{ni)^<wlGrRe + 0{T-^) . (1.8) 

With some very mild technical restrictions on the forcing, Doering and 
Foias [19j then showed that Navier-Stokes solutions obey Gr <cRe^ . This 
turns (jl.Sp into 

{Qi)^ < cTulRe^ + O (T-i) . (1.9) 

In fact is the time-averaged energy dissipation rate per unit vol- 

ume over [0, T] and allows us to form and bound from above the inverse 
Kolmogorov length scale A^^ 

A-4 = !:i^ ^ LA,i<ci?e3/4 + 0(T-i/^). (1.10) 
An estimate for the inverse Taylor micro-scale A^,^^ can also be found from 

(lESl) 

|2\ \l/2 

<ci?ei/2 + o(r-V2). (1.11) 

Both these upper bounds gratifyingly coincide with the results of statistical 
turbulence theory [U El [3] although the fact that they are bounds allows 
for structures to occur in a flow whose natural scales are larger [20j. The 
question is now clear : can we construct and hound from above a sequence of 
inverse length scales associated with the higher moments ? 




^Doering and Foias [19] took narrow-band forcing around a specific wave-number but 
a wave- number spectrum which is cut off both above and below is sufficient. 
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2 A scaling property and length scale estimates 

Leray's energy inequality (jl.7p is valid for weak solutions and thus the es- 
timate ()1.9p is rigorous, although the existence and uniqueness of solutions 
for arbitrarily long times remain an open problem. While it is possible to 
subscribe to the view that difficulties in flow resolution could be a symptom 
of the lack of uniqueness of weak solutions, in tandem it ought also to be 
acknowledged that these difficulties may simply be caused by the practical 
challenges of working on a system where even the naturally largest scale 
(other than L) lies close to the limit of what can currently be resolved. The 
spirit of this paper is such that results on weak solutions are assumed to 
be sufficiently physical to reflect the reality of turbulent flows, provided T 
is taken large enougtEI. This strategy allows the estimation of an infinite 
hierarchy of time averages of powers of the for weak solutions on [0, T] 
without having to appeal to point-wise estimates that the solution of the 
regularity problem would require. In turn, these time averages allow us to 
define and explore the natural length scales inherent in the system. The fol- 
lowing result was proved in |21] under the assumption that strong solutions 
exist. Here it is demonstrated for weak solutions : 

Theorem 2.1 Weak solutions of the three dimensional Navier-Stokes equa- 
tions satisfy 

>(r^f^m)°'")^<ci?e3 + 0(T-^), l<m<oo, (2.12) 

where c is a uniform constant and 

'2i7n 

Remark : The exponent Q!m = 2?7i/(4m — 3) within ()2.12p appears as a nat- 
ural scaling of the Navier-Stokes equations, consistent with the application 
of Holder and Sobolev inequalities. Note that when m = l the value ai = 2 
is consistent with (11.91). 



Proof : The proof is based on a result of Foias, Guillope and Temam [22j 
(Theorem 3.1) for weak solutions which, when modified in the manner of 
Doering and Foias [19] . furnishes us with the following time averaged esti- 
mate 



H^""-') <CNL-^u^^Re^ + 0{T-^) , (2.14) 

T 



^While existence and uniqueness of solutions is easily proved for small values of T 
[15| 1161 117| , larger values than this are necessary to make sense of long-time averages. 
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where 



Hn= [ \V^ufdV= [ k^^\u\^Sk, (2.15) 

where Hi = J^lS/ul^ dV = J^IljI^ dV . An interpolation between ||a^|[2m and 
|[a;||2 is found using Hn 

M2m<CN,m\\V''-^ur2\Ml-\ "= 2m(Ar~-l) ' ^^-^^^ 

for > 3. ||w||2m is now raised to the power Om, which is to be determined, 
and the time average of this is estimated as 

/IL ill^^iN ^ / II 11/ 1l|(-'-~'')"™■ 



T 



2N-1 \ U2^ a;am 



/ T 

1 \ ^aam{2N-l) i (i-a)c.„ \ \-\aa^{2N -\) 

< ^%^m{HN-' ) ™{2iV-i) \ (2.17) 



T 

An exphcit upper bound in terms of Re is available only if the exponent of 
Hi within the average is unity ; that is 

(l — a)am 

^ ' -1. (2.18) 



2-aQ„(2A^-l) 
This determines a^, uniformly in A^, as 

2m 

a„ = j— (2.19) 

Using the estimate in (j2.14p . and that for {Hi), the result follows. The 
constant CN,m can be minimized by choosing A^ = 3. C3,m. does not blow up 
even when m = oo; thus we take the largest value of Cg ^ and call this c. ■ 

3 Definition of the inverse length scales 

Based on the definition of the inverse Kolmogorov length A^^ in (jl.lOp a 
generalization of this to a hierarchy of inverse lengths suggests the 
definition 

[LX-^^f--:={{w,'n^r-)^. (3.20) 
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The Xm are interpreted as resolution lengths in the space-time averaged 
sense for 1 < m < oo : 

Many turbulent structures have natural inverse gross length scales lying in 
the range between i?e^/^ and Re^^^, but crinkles forming at finer scales may 
ultimately grow to be dominant and then become the cause of resolution dif- 
ficulties [HElEKiaill]. For m>l the A m are interpreted here as the length 
scales corresponding to these deeper intermittent events. The upper bounds 
displayed in (|3.2ip . as the Table shows, range from the inverse Kolmogorov 
length Re^^^ at m = l to Re^ for m = oo. Computationally it may be hard 
to get far beyond m = l: for example, m = 9/8 corresponds to Re^, which 
is close to modern resolutions even for modest values of Re. Thereafter the 
rise in the exponent 3/2am is steep. Indeed, in the very high m limit, the 
i?e^ bound has an exponent four times greater than the Kolmogorov length ; 
this lies well below molecular scales where the equations are invalid. 



m 


1 


9/8 


3/2 


2 


3 




oo 


3 /2am 


3/4 


1 


3/2 


15/8 


9/4 




3 


dm 


3 


2 


1 


3/5 


1/2 








Table 1: Values of the Re-exponent 3/2q;„j = 3 (l — ^) , and dm = 



An interesting question is how the existence of this continuum of finer scales 
might be interpreted physically? To do so rigorously without a regularity 
proof is difficult but a very informal physical interpretation is possible in 
terms of the familiar concept of a cascade. One of the simplest cascade 
models is the so-called /3-model of Prisch, Sulem and Nelkin |23| who, in 
analogy with Mandelbrot's ideas [23], modelled a Richardson cascade by 
taking a vortex of scale io^L and then allowed a cascade of daughter vor- 
tices, each of scale in- The idea was based on domain halving at each step 
such that io/in = 2^. The self-similarity dimension d (similar to fractal di- 
mension) was then introduced by considering the number of offspring at 
each step as 2^^ : 2 for the halving of a line ; 4 for the halving of each direc- 
tion in the plane ; and likewise 8 for the cube, d is then formally allowed 
to take non-integer values. In d dimensions the corrections to the usual 
Kolmogorov scaling calculations for velocity, turn-over time etc appeared as 
multiplicative factors proportional to {£o/in)^^~'^''^^ '■ see [ll[23]. Equating 
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the turn-over and viscous times in the standard manner one arrives at {id 
is their viscous dissipation lengthj^ 

io/id^Re^. (3.22) 

This gives the usual Kolmogorov inverse scale of i?e^/^ in a fully three- 
dimensional domain but is shifted upwards for smaller values of d. Taking 
this idea as our physical analogy we compare (j3.22p to the upper bound in 
([3:2T]1 to get 

3 

dm + l = 2am =^ dm = - r, (3.23) 

4m — 3 

where an m-label has been appended to d. Thus we are able to assign 
a corresponding self-similarity dimension d^ to lower-dimensional vortical 
structures that require values of ?7i > 1 for their resolution. Note that d-m 
never goes negative. Models more sophisticated than the /3-model, such 
as the bi- fractal and multi-fractal models [Ij, are more difficult to use as 
analogies as they would require data fitting. 
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